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^  A  common  problem  in  controlling  a  large  flexible 
structure  with  a  reduced  order  model  is  the  excitation  of 
the  residual  (non-model led)  mo^s  by  the  control  inputs 
resulting  in  control  spillove:^  In  this  thesis  a  new 
approach  to  eliminate  controy spillover  is  examined. 

This  approach,  based  on  the  independent  modal-space  control 
method,  uses  a  finite  number  of  spatially  distributed  input 
points  to  eliminate  control  spillover.  Control  of  an 
undamped  beam  and  a  finite  degree  of  freedom  truss  is 
accomplished  through  con^uter  simulation.  The  results  of 
the  seven  example  problems  shows  that  this  new  approach 


does  effectively  eliminate  control  spillover. 
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INTRODUCTION 


With  the  Increased  capahllltles  of  the  space  shuttle 
man  will  now  he  able  to  build  large  structures  In  space. 
These  structxires  will  have  applications  In  many  areas  that 
will  demand  precise  pointing  and  control  such  as  solar 
power  stations,  radio  astronomy,  and  comminl  cat  Ions 
antennas.  In  order  to  satisfy  mission  performance  require¬ 
ments,  these  large  structxires  will  require  control  systems 
that  will  actively  suppress  the  vibrations  of  the  structure 
(IJ. 

As  a  space  structure  Increases  in  size  the  complexity 
of  the  control  problem  Is  Increased  by  the  number  of  modes 
that  must  be  controlled.  However,  due  to  many  limitations 
it  is  not  always  possible  to  control  all  of  the  modes  of 
the  structure  and  a  reduced  order  model  of  the  flexible 
system  must  be  used.  But  this  reduced  order  model  creates 
new  problems  since  the  residual  (non-modeled)  modes  of  the 
structure  can  be  excited  by  the  control  Inputs  causing 
control  spillover.  The  reduction  and  elimination  of  this 
spillover  is  of  critical  importance  In  designing  a  precise 


control  system. 


In  this  paper  we  will  examine  a  new  approach  forwarded 
in  [2]  that  is  based  on  the  independent .modal -space  control 
method.  This  approach  uses  a  finite  number  of  spatially 
distributed  input  points  that  essentially  eliminates 
control  spillover  on  the  modes  beyond  the  control  model. 

We  will  first  give  an  overview  of  the  formulation  of  this 
method,  then  we  will  apply  the  approach  through  computer 
simulation  to  an  undamped  beam  and  an  undamped  truss 
problem.  We  will  use  the  results  to  determine  if  the 
approach  to  the  control  problem  proposed  in  [2]  cam  reduce 
or  eliminate  the  control  spillover. 


CHAPTER  I 


COHTROL  SYSTEM  FORHULATIOH 


Introduction 

In  tills  chapter  we  will  overview  the  formulation  for 
the  control  method  presented  in  [2].  The  flexible  struc¬ 
ture  will  be  examined  as  a  distributed-parameter  system. 

By  discretizing  this  system  in  space  one  converts  the 
partial  differential  equation  into  an  infinite  set  of 
ordinary  differential  equations  for  the  modal  coordinates. 
By  using  the  independent  modal -space  control  method  (IHSC) 
and  the  integral  formulation  of  the  structural  dynamic 
system  one  can  design  a  desired  control  system. 

The  Distributed-Parameter  System 

The  displacement,  u(x,  t),  of  an  arbitrary  point  x  of 
an  undamped  nongyroscopic  distributed- parameter  system  can 
be  expressed  by  a  partial  differential  equation  [3] 

Lu(x,  t)+M(x)d2u(x,  t)/dt2  =  f (X, t)  (1 

which  must  be  satisfied  at  every  point  x  of  the  domain  D. 
In  equation  (1)  L  is  a  self-adjoint,  positive  definite 


linear  differential  operator  expressing  the  system  stiff¬ 
ness,  H(x)  is  the  mass  distribution,  and  f(x, t)  is  a 
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distributed  input.  In  addition  the  displacement  must  also 
satisfy  given  boundary  conditions. 

When  solving  equation  (1)  it  is  observed  that  the 
system  possesses  a  denumerable  Infinite  set  of  eigenvalues 
Xp,  ¥rhich  are  related  to  the  natural  frequencies  by 

r=l,2, ....  and  eigenfunctions  ^p(x).  Because  L  is 
self-adjoint  the  eigenfunctions  are  orthogonal  and  can  be 
normalized  to  satisfy 


I 


p  db=<5  . 

j)  s  r  rs 


r,  s:  1.2,.. 


p  L.p  dD=u  ^6  I  r,  s:l,2,... 
jj  s  r  r  rs 


(3) 


(3) 


in  which  <Sp5  is  the  Kronecker  delta. 

By  using  the  expansion  theorem  [4]  the  solution  for 
the  system  can  be  represented  by  an  infinite  series  of 
space -dependent  eigenfunctions  ^p(x)  multiplied  by  a  time- 
dependent  generalized  modal  coordinate  qp(t)  of  the  form 
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u(x,  t)=  2  ^r(*)<Ir<t). 
r=  1 


(4) 


By  placing  equation  (4)  into  equation  (1),  multiplying  by 
^s(x),  Integrating  over  D,  and  using  the  orthogonality 
equations  (2)  and  (3)  we  obtain 

qp(t)+Wp2q^(t)=fp(t),  r=l,2,...  (5) 

where 


f  (t)=  p  (x)f  (x,  t)dD, 
r  n  V 


m 


i 


§ 


V 

»  *Jm  ‘ 


SI 


V-V. 

is 


mi 


r=l,2,.  .  . 


(6) 
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are  generalized  modal  control  forces.  Because  It  is 
seldom  feasible  to  control  ttie  entire  Infinity  of  modes  we 
must  reduce  the  distributed-parameter  model  to  a  finite 
number  of  modes  that  are  sufficient  to  accurately  model 
the  system. 

Independent  Modal -Space  Control 

If  the  mathematical  model  of  our  system  Is  very  large 
It  Is  not  always  practical  and  economically  feasible  to 
control  all  of  the  modelled  modes.  Therefore,  the  model 
Is  broken  Into  two  groups  the  controlled  modes  and  the 
uncontrolled  or  residual  modes.  By  assiimlng  that  the 
control  Input  Is  Implemented  by  m  point  actuators,  then 
the  actuator  forces  can  be  described  as  discrete  by  [4] 

m 

f  (X,  t)=  Z  Fi(t)<J(x-Xi),  (7) 

1=1 

where  dtx-x^)  Is  a  spatial  Dirac  delta  function.  By 
placing  equation  (7)  into  equation  (6)  we  obtain 

m 

fr(t)=  S  ^r(»l)Fl(t),  r=l,e .  (8) 

1=1 

By  examining  equation  (6)  one  can  see  that  the  actual 
control  forces  F^(t)  can  excite  each  mode.  However,  these 
forces  are  designed  to  suppress  vibration  in  the  controlled 
modes,  which  results  In  excitation  of  the  uncontrolled 
modes.  'Rils  excitation  of  the  uncontrolled  modes  Is 
referred  to  as  control  spillover.  This  spillover  exists 


V  *..  •*« 


because  the  actuators  are  discrete  elements.  If  distri¬ 


buted  actuators  were  available,  the  modal  control  forces 
could  be  generated  so  that  they  do  not  excite  the  residual 
modes  and  therefore  eliminate  control  spillover. 

We  can  put  equations  (5)  and  (d)  into  compact  form  by 


defining  the  following  vectors 

q(t)  =  (qi(t).  qgCt).  .  .  .  jT,  (9: 

f(t)  =  {fi(t).  f2(t),  .  .  .  )T,  (lO; 

F(t)  =  (Fi(t),  F2(t) . Fni(t)jT  (li; 

and  the  matrices 

O^sdiagCwj^,  ...),  (12! 

Bs  (Xs)l(  rsl,2,....  S:l.&|...m.  (13 

Then  equation  (5)  can  be  expressed  as 

q(t)+Q2q(t)rf(t)  (14 

and  equation  (6)  becomes 

f(t)=BF(t).  (15 


From  equation  (15)  it  can  be  seen  that  the  actual 
control  force  can  be  obtained  from  the  modal  control  force 
by 

F(t)  =B-*f (t).  (16 

The  modal  control  forces  in  general  can  be  designed  as 
desired  in  the  form  of  fx>(qiq)  to  design  a  closed-loop 
control  system.  However,  if  each  modal  force  fp  is  a 
feedback  function  of  all  modal  coordinates  q,  then  the 
modal  responses  become  coupled  during  control,  so  that 
modal  coordinates  q^  no  longer  represent  natural 


coordinates  for  the  controlled  system.  On  the  other  hand, 

If  each  modal  Input  is  designed  as  a  function  In  the  form 
of  feedback  of  only  Its  own  modal  coordinates,  1.  e. 
fr(qr>^r)  then  the  response  of  each  controlled  mode  will  be 
Independent  of  the  responses  of  other  modes  and  each  qp 
will  retain  its  natural  property.  This  particular  tech¬ 
nique  of  modal  feedback  control  is  known  as  Independent 
Modal -Space  Control  [5]. 

Integral  Formulation  of  Structural  Dynamic  Systems 

The  structural  dynamic  equation  of  a  self-adjoint, 
positive  definite  distributed-parameter  system  can  also  be 
written  alternately  in  the  integral  form  [6] 

u(x,  t)  =u(x,  t)  ♦|m(<n)C(x, -nlucn.  t)dii=|c(x,  Ii)f  (t),  t)dii,  (17) 

where  x,  t)  denote  the  structural  domain.  u(x,  t)  is  the 
instantaneous  static  deflection  caused  by  f('ii,  t).  C(x,  ii) 
is  the  flexibility  Influence  function  that  is  dependent  on 
the  boundary  conditions,  material  properties,  and  geometry 
of  the  system. 

It  can  be  shown  [2]  that  if  we  express  the  spatial 
input  distribution  f(‘n,  t)  and  the  flexibility  function 
C(x,  Ti)  in  the  form 

n 

f  Cn,  t)=  2  mCn)^p(n)fp(t), 
r=  1 


(18) 


fi 


C(x,  T1)=  S 
r:  1 

where  are  eigenfunctions  and  fp(t)  are  modal  input 
forces,  then  the  Instantaneous  static  displacement  will 
have  the  form 


u(x, t):  2  ^p(x)«r-efp(t) 
r:  1 

If  we  introduce  the  vector 


r=  1, 2, .  .  .  n. 


...  ^r(x)  •••  ^n<x)J 
then  equation  (20)  can  be  written  as 
u(x,  t)=bHTQH-2f (t). 


where 


f (t)  =  (fi  #2  .  .  .  fr  •  •  •  fnJ- 


If  we  now  look  at  equation  (17)  and  consider  u(x, t) 
at  a  finite  n\imber  of  points  with  the  load  f(x,  t)  at  a 
finite  number  of  input  points  and  using  the  collocation 
technique  [2]  so  that  the  controlled  points  are  the  same 
as  the’  input  locations,  then  we  can  rewrite  equation 

(17)  as 


ic  c  ...c  ;c  ...c 

11  12  Im  •  Im'fl  loo 


c  ;  c  .  .  .  c 
.  mm_  _  ann±i_  _  _  a>» 


m+ii 


In  equation  (24) 


Cij  =  C(Xi,  xj) 


1.  J  =  l.  2, 


V 


•  J 
*  .  .‘.‘1 
V* 


Ui  =  u{Xi,  t),  i=l,2.  .. 

f  (Xk.  t)  :Fn(t)<5(x-Xit)  K=l,2,...in 

and  It  can  be  easily  rearranged  as 


u  " 

1 

-  • 

c  c  .  .  .  c 

11  12  im 

u  * 
_m_ 

c  .  c 

-ml  _____  aim_ 

u 

m-f  1 

c  . 

m-t-l  1 

c  . 

.  ®1  J 

(26) 

(27) 


(26) 


where  »  denotes  a  desired  (prescribed)  behavior.  Hence 
u”  (x^^)  l=l,2....m,  Is  the  desired  Instantaneous  static 

displacement  for  the  controlled  points  Xj^.  1  =  1,2,. ..m, 
which  can  be  realized  by  designing  the  Input  forces  F.  The 
C  matrix  contains  all  the  flexibility  coefficients  for  the 
controlled  points  and  the  Cp  matrix  contains  the  coeffi¬ 
cients  for  the  residual  points.  From  equation  (26)  we  can 
easily  find  the  unique  solution 

F(t)=C‘lu",  (29) 

where 

u^iCui*  ug"  ...  ^n,*)'^.  (30) 


Specified  Response  Model 

In  view  of  equation  (29)  we  need  to  devise  a  method 


by  which  the  desired  Instantaneous  static  response  vector 
u"  can  be  specified.  From  equation  (4)  it  is  easy  to  show 
that  a  desired  displacement  profile  u*(x,  t)  can  be  express- 


ed  in  a  similar  form,  but  by  using  a  desired  modal  coordi¬ 
nate  q” 

n 

u“(x.  t):  S  ^r(x)qr*‘(t).  (31) 

r=l 

It  tben  follows  that  we  can  state  a  set  of  second  order 
model  dynamics  for  the  coordinates  q^"  as 

qr"+Qr^qr*  =  fr"(^)*  r:l,a,...n,  (32) 

vdiere  fp”  is  a  model  forcing  function.  It  cam  be  seen  that 
by  adjusting  the  model  coordinates  q^*  in  eqniation  (31)  one 
can  generate  a  desired  displacement  u”(x,  t).  This  can  be 
accomplisbed  by  designing  tbe  model  inputs  in  the  form 
fr*  (‘Ir"* ‘Ir" )  sucb  tbat  qp"  bebaves  as  desired. 

By  using  the  same  method  as  used  in  producing  equation 
(22)  it  can  now  be  shown  that 

uMx,  t)=bHT0H‘2fN"(t).  (33) 

If  we  now  evaluate  equation  (33)  at  a  finite  number  of  m 


points,  u*‘(Xi):Uj^*  l:l,2,...m,  and  form  the  vector  u"  we 
obtain 

u"=  BuTnjj-2f«  (t)  (34) 

®N=C^(*l)  )^(*2)  •••  ^M(*m))  (3®) 

bN(XK)=  (^i(Xk)  ^2(*k)  •••  ^r(*K)  •••^n(*K)j'^  (36) 

If  we  now  substitute  equation  (34)  into  equation  (29)  we 
obtain 

F(t)=C-*u"=C-lBN'^%-2f„»(t)  (37) 


Vhere  f|i"(t)  is  an  n-vector  of  model  inputs  fii*(t)  =  {fi* 
f2*  ...  designed  by  independent  modal-space  control 


11 


m 


of  the  model  given  by  equation  (32).  This  is  equivalent  to 
designing  a  model  following  input  force  F  for  the  actual 
structure  in  an  open-loop  sense  [2]. 

Open  Loop  Control 

If  the  input  forces  F  in  equation  (15)  are  designed  on 
the  basis  of  model  inputs  f*»  equation  (37),  of  a  (modal) 


reference  model,  equation  (32),  and  are  applied  directly  to 
the  structure  through  equation  (14)  without  regard  to  the 
actual  response  of  the  structure,  then  the  response  design 
is  essentially  an  open-loop  control  system.  By  introducing 
equation  (37)  into  equation  (14),  which  is  the  structural 


,;vv 


modal 

set,  we  obtain 

iV! 

q ( t ) +o2q  ( t )  =  BgC- 1 BN'>^QN"2f n"  ( t ) , 

(30) 

where 

®s=  C®n'^  1 

(39) 

BR=tbR(Xi)  ...  bR(Xni)], 

(40) 

»>R(Xk)  =  t^n+l  (*k)  ^n+2<Xk) 

(41) 

*'  / 

and 

qMqN*^ 

(42) 

c 

qN={qi  q2  •  •  •  qnl"^*  ‘lR=l‘In+l  ‘In+2  •••)'*’• 

(43) 

>v! 

If  we 

expand  equation  (38)  in  terms  of  Br  and  Bp  we 

q(t)+02q(t)  =  [EpT  BpTjTc-lBpTo^-SfHMt). 

get 

(44) 

One  can  easily  see  that  there  will  be  control  spillover  to 
the  residual  dynamics  qH(t)  from  the  model  input  f*. 
However,  if  the  structure  had  a  finite  number,  N,  of  modes 


and  63:611  then  there  would  be  no  residual  modes  to  incur 
spillover. 

Another  approach  to  eliminate  spillover  to  residual 
modes  r^n-fl,  n-^S, .  .  .  if  the  system  has  an  infinite  number 
of  modes,  is  to  increase  the  number  of  inputs  m— »cd.  It 

can  be  shown  [7,21  that  the  triple  matrix  product  BnC”^6||7' 
will  diagonalize  to  the  eigenvalues  matrix  Oji^.  This  can 
be  seen  by  remembering  that  the  eigenvalue  problem  of  a 
self-adjoint  distributed- parameter  system  has  the  orthogo¬ 
nality  property 

(#  (x)L#  (x)dD=«  £<5  ,  (45 

s  r  r  rs 

where  L  is  the  distributed  stiffness  operator.  6y  recall¬ 
ing  the  form  of  6k  and  Bp,  noting  that  the  C~^  matrix  is 
the  stiffness  matrix  1C  and  considering  the  distributed 
orthogonality,  one  can  see  that  the  triple  matrix  product 
BkC'^Bk*^  is  the  discrete  approximation  of  equation  (45). 
Therefore  as  the  nximber  of  inputs  goes  to  infinity  the 
system  starts  to  resemble  the  distributed-parameter  system 
more  closely  which  results  in 

BHC‘‘BH'r=6HKBM'f  - >  [Qh^I  ("^6 

BrC‘^BhT=6rK6hT  - »  0.  (47 

If  we  apply  this  to  equation  (44)  we  get 

q(t)+Q2q(t)=[I50lTf«(t)  (48 

where  [I]  is  the  NxN  Identity  matrix.  This  is  essentially 
the  same  as  equation  (14).  Since  the  modal  inputs  f”(t) 


are  designed  in  the  form  f^Cq^iq”).  provided  that  the 
Initial  conditions  on  the  model  and  the  structure  are  the 
same,  the  open-loop  approach  for  modal  model  following  Is 
equivalent  to  a  modal  feedback  control  law  on  the  actual 
structure  since  one  will  have  qrq*  tlO. 

Closed-Loop  Control 

Por  the  closed-loop  case  the  Input  forces  t* (t) 


are  designed  In  the  form  f*(qHiqN) 
fll" 

where  and  Hg  are  HxH  modal  control  gains  of  the  font 

Hi=PtIl.  H2=a[I].  (50 

If  we  apply  this  to  the  model  n-modes  of  equation  (38)  we 
get  the  form 

We  can  rearrange  this  equation  and  we  get  the  form 

qiI-aBnC“  ^Bjj;'^Oij"2qjj+  (Oj|2-pBj|C"  ^Bjj'^Djj"^)  qN=0.  (52 

It  can  he  easily  seen  that  the  residual  dynamics  will  have 
the  form 


If  we  examine  eqpiatlons  (52)  and  (53)  as  the  number  of 
inputs  Is  Increased  to  Infinity  It  can  easily  be  seen,  by 
virtue  of  equations  (46)  and  (47),  that  the  control  model, 
equation  (52),  goes  to  the  Ideal  case  where  the  controlled 


distributed-parameter  system  behaves  as 


(54 


and  the  residual  excitation  totally  disappears.  Essential¬ 
ly,  the  controlled  behavior  exhibited  by  equation  (54)  is 
precisely  the  behavior  that  will  be  obtained  if  a  spatially 
continuous  input  as  in  equation  (18)  is  used  with 

fp=aqp4^Eqp  r:l,£,  ...n.  (55) 

By  choosing  a  as  a  real  negative  number  and  B  as  a  real 
number  less  than  equation  (54)  represents  a  stable 

system  where  a  can  be  regarded  as  a  modal  damping  coeffi¬ 
cient  and  B  as  a  modal  stiffness  coefficient.  With  the 
elimination  of  the  residual  modes  the  n  controlled  modes 
will  be  uniformly  damped  and  stiffened. 

We  have  seen  that  the  feedback  control  of  the  struc¬ 
ture  is  performed  via  an  input  control  force  in  the  form 

F(t)=C-iBN'^OH-2fH(t).  (56) 

It  can  be  shown  [2)  that  the  spatially  discrete  inputs  F 
are  redistributed  by  synthesizing  the  modal  feedback 
quantity  fH<t)  through  the  C  matrix  so  as  to  create 
effective  continuously  distributed  damping  and  stiffness 
Inputs  (see  [2]  for  proof).  Therefore  equation  (56)  is 
referred  to  as  "spatial  modal  input-distribution  control" 
[2]  that  has  an  inherent  spatial  orthogonal  non-pass 
filtering  property  with  respect  to  the  residual  modes. 

This  occurs  through  the  gain  factor  that  eliminates 

spillover  as  the  n\imber  of  inputs  increase.  However,  the 
number  of  inputs  for  most  structures  does  not  have  to  go  to 
infinity  before  the  model  essentially  reaches  the  ideal 


CHAPTER  II 


CONTROL  OF  AN  UNDAMPED  BEAM 

Introduction 

In  this  chapter  the  closed-loop  feedhacK  control 
technique  developed  In  the  last  chapter  will  he  applied  to 
a  slnqply  supported  beam.  The  method  of  linplementlng  the 
control  system  will  be  examined,  then  three  examples  of 
the  control  system  will  be  shown. 

Implementation 

To  Implement  the  closed- loop  feedback  control  method 
we  must  put  the  modal  ecjuatlons  Into  state  form.  Ihe 
modal  equation  for  a  distributed-parameter  system,  equation 
(5),  can  easily  be  put  Into  state  form  by  using  the 
auxiliary  variable  lip(t)  defined  by 

qr(t) =«rtr(^) •  r:l,2 .  (57) 

The  resulting  equation  In  modal  state  form  Is  found  to  be 

[ai 

Wi.(t)=ArWi,(t)*Wr(t),  (58) 

where  Wp(t)  is  the  modal  state  vector,  Wp(t)  the  associated 
modal  control  vector,  and  Ap  the  coefficient  matrix  defined 


Wr(t)  =  [qj.(t) 

Wj,(t)  =  [0  fp(t)/WrlT, 


0  w 

A  =  ^ 

r  .w  0 

r 


r=l,  2, . . . 


However,  the  Independent  modal-space  control  method 
designs  the  control  vector  in  the  form 

Wp(t)=QpWp(t),  (62) 

where  Qp  is  a  2x2  modal  control  gain  matrix.  nie  ideal 
modal  gain  can  be  found  by  putting  equation  (54)  into  state 
form  since  it  represents  the  distributed-parameter  system 
behavior  for  the  ideal  control  model.  nie  resulting 
control  design  model  has  the  form 

Wp(t)  :ApWp(t)  ♦GpWp(t)  t  rsl,2,  ...n  (63) 

where  the  ideal  gain  matrix  is  defined  by 


O  0 

G  : 

p/w  a 


r=  1,  2, .  .  .  n. 


The  ideal  closed-loop  eigenvalues  can  be  found  from 
equation  (52)  and  will  have  the  form 

ip=<x/2  i  via2-4«p2-4p/2.  (65) 

To  implement  the  control  design  the  continuous -time 
modal  state  equations  must  be  discretized  by  sampling  in 
time  to  permit  digital  state  estimation  and  control.  This 
can  be  shown  [9]  to  result  in  the  form 


Wp(k-^1 )  :$pWp(k)  ■frpWp(k).  ksl.2,...,  rsl,2.  ..n  (66) 


where 


5^(r)  = 


r  (T)=(i/w  ) 
r  r 


cosfa)  T  sinu  T 
r  r 

-Sln(i)  T  cosw  T  , 

L  r  r  J 

[sinu  T  1-cosu)  t’ 

r  r 

-(1-cosw  T)  sinw  T  , 

r  r  J 


in  whicli  T  denotes  tlie  sampling  time  and  the  arguments  k^l 
and  K  denote  the  times  t=(K4^1)r  and  t:kT  respectively. 

The  simply  supported  beam  (  Fig.  1  )  is  chosen  with 
tinit  length,  unit  mass  distribution,  and  unit  stiffness. 
The  eigenfunctions  and  natural  frequencies  of  the  beam  are 
known  to  be  [10] 


j2/inL  sinnrx/L, 
Wi»sr2w2  JeI/boL^, 


rs 1,  2, . . . 


where  m  is  the  mass,  L  is  the  length,  and  El  is  the 
stiffness.  The  closed  form  of  the  influence  function 
C(x,  t))  for  the  beam  with  E  inputs  can  be  fotind  to  be  [2,9] 
C(Xi,  Tii)  =Xi/6EI  ( Xi2,,^/L+<Hi3/L+2'»iiL-Xi2- 31,^2 ) ,  Xi‘T|i  (71 
for  i=l,2,...E.  If  then  the  ii  and  x  termis  are 

interchanged. 


f(nk.t) 


Fig.  1.  Simply  supported  uniform  beam 
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Computer  simulation  of  the  control  system  is  accom¬ 
plished  by  the  use  of  the  coniputer  program  PHTCHT  (Appendix 
A).  The  program  output  consists  of  plots  of  the  displace¬ 
ment  response  of  a  chosen  control  point  versus  time  for  the 
reduced  order  model  and  the  structxire  which  Includes  the 
residual  modes.  The  input  force  at  the  same  point  versus 
time  is  also  plotted.  nie  horizontal  axis  on  all  plots 
reflect  time  and  the  vertical  axis  reflect  either  displace¬ 
ment  or  the  input  force.  Since  a  beam  has  an  infinite 
number  of  modes  the  structure  is  modelled  by  a  finite 
nximber,  N^,  of  modes  that  are  sufficient  to  accurately 
represent  the  Infinite  set  of  modes.  The  reduced  order 
model  that  is  controlled  consists  of  Hsn  modes  where  H<Hs 
and  the  number  of  inputs  will  be  referred  to  as  E.  The 
input  points  are  picked  arbitrarily  along  the  beam,  but  no 
points  are  allowed  to  coincide  with  any  of  the  structures 
nodes  since  that  would  reduce  their  effectiveness.  The 
initial  modal  state  disturbances  for  all  the  examples  are 
chosen  as  qp(0)r0  and  4^(0):  1.0  units/second,  r=l,2,...N. 

Example  1 

In  this  first  example  the  beam  has  unit  length,  L=l, 
and  the  structure  is  represented  by  ten  modes,  Ns=10.  The 
reduced  order  model  consists  of  the  first  four  modes  of 
the  structure,  N=4.  The  damping  coefficient  is  chosen  as 
a: -40  and  the  stiffness  coefficient  is  chosen  as  p=0.  The 


£0 


point  to  observe  tbe  model  and  structure  responses  at  was 
chosen  as  x:0.  21  units.  The  time  tics  on  the  plots  are  at 


0.  5  second  Intervals  and  the  displacement  and  force  tics 
are  at  1.0  unit  Intervals.  The  displacement  response 
plots  for  both  the  model  and  the  structure  can  be  seen  in 
Figures  2-11.  Vhen  £=2  the  input  points  are  at  x=0.  09  and 
z=0.  21.  When  K=4  the  points  x=0.  31»  0.42  are  added  as 
input  points.  When  C=6  the  points  x=0.  53,  O.  64  are  added, 
for  i:=6  points  x=0.  73,  0.64  are  added  and  for  K=9  the  point 
x=0.  93  is  added  as  an  input  point. 

As  we  examine  Plgtires  2-11  it  can  be  seen  that  as  the 
number  of  inputs  Increase  the  model  and  the  structure 
responses  start  looking  more  and  more  alike  until  they 


match  each  other  when  K:9.  If  we  compare  Figures  10  and 
11  to  Figure  12,  which  is  the  response  of  the  ideal 
model,  equation  (54),  that  is  based  on  an  infinite  number 
of  input  points,  1.  e.  on  a  continuously  distributed  input  we 
see  that  they  are  identical.  We  can  confirm  that  the 


actual  model  approaches  the  ideal  model  by  comparing  the 
gain  matrices  and  the  closed- loop  eigenvalues  of  the 


actual  model  to  the  ideal  gain  matrix  and  eigenvalues.  The 
ideal  gain  matrix  has  the  form 


G: 


(72) 


When  K=2  the  non-ideal  gain  matrix  is 


0 

0 

0 

0 

0 

0 

0 

0 

0  - 

34.  37 

0 

-12.  14 

0 

-3.  46 

0 

.  91 

0 

0 

0 

0 

0 

0 

0 

0 

0  - 

12.  14 

0 

-7.  03 

0 

-5.  49 

0 

-4.  62 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-3.  48 

0 

-5.  49 

0 

-6.  98 

0 

-7.  60 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  91 

0 

-4.  62 

0 

-7.  60 

0 

-8.  95 

the 

non-: 

Ideal 

gain  matrix 

is 

0 

0 

0 

0 

0 

0 

0 

0 

0  - 

39.  92 

0 

-.  45 

0 

-.  93 

0 

-.  97 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.  45 

0 

-37.  22 

0 

-5.  72 

0 

6.  13 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.  93 

0 

-5.  72 

0 

-27.  82 

0 

-13.  48 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-.  97 

0 

6.  13 

0 

-13.  48 

0 

-23.  11 

'  the 

non-: 

ideal 

gain  matrix 

Is 

0 

0 

0 

0 

0 

0 

0 

0 

0  - 

40.  00 

0 

-.  0002 

0 

.  0004 

0 

.  002 

0 

0 

0 

0 

0 

0 

0 

0 

0  - 

.  0002 

0 

-40.  00 

0 

.  004 

0 

.  0001 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  0004 

0 

.  004 

0 

-39.  92 

0 

.  03 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.  002 

0 

,  0001 

0 

.  03 

0 

-39.  66 

Table  1.  Example  i  closed-loop  eigenvalues 


IDEAL 


-2.  60 
-37.  40 
■20.  0ii34.  04 
•20.  0!i66.  55 


-3.  08 
-33.  56 
-2.  5ii38.  35 
-3.  41168.  65 


K=6 


-2.  61 
-37.  32 
16.  71i34.  89 
13.  91i86.  00 


K:9 


-2.  60 
-37.  39 

-20.  01134.  04 
-20.  01186.  55 


-20.  011156.  6 


-4.  511157.  5 


-11.  411156.  5 


19.  811156.  7 


By  examining  tb.ese  matrices  and  Table  1  It  becomes  obvious 
that  as  the  nuinber  of  Inputs  Increase  the  actual  model 
does  approach  the  Ideal  model  and  the  effect  of  the 
residual  dynamics  on  the  structural  response  Is  practically 
eliminated.  Another  Item  to  note  Is  the  Input  force.  In 
Figures  14-17  It  can  be  seen  that  as  the  number  of  Inputs 
Increase  the  Input  force  decreases.  This  Is  logical  since 
the  amoxmt  of  force  required  to  control  the  system  Is  being 
distributed  to  more  Input  points  and  the  force  response 
eventually  matches  the  Ideal,  Figure  13. 

One  more  point  of  Interest  would  be  to  see  the  effect 
on  the  model  and  the  structure  response  for  K:9  If  the 
number  of  sufficient  modes  to  represent  the  structure  Is 
Increased.  In  Figures  18  and  19  the  structure  l.s  in¬ 
creased  from  Ns=10  to  KgslS  and  the  effect  Is  negligible. 


Fig.  2.  Example  i  model  dlsplacement/tlme  K:Z 


Fig.  3.  Example  1  structure  displacement/time 
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Example  2 

llils  closed-loop  example  is  very  similar  to  Exanqple  1 
since  L=l,  NgslOi  E=0,  and  the  observed  point  will  remain 
as  x=0.  21.  However,  we  will  increase  the  reduced  order 
model  to  N=6  and  the  damping  coefficient  is  increased  to 
a: -80.  The  time  tics  and  the  displacement  and  force  tics 
will  remain  the  same  as  in  Example  1.  All  the  Input  points 
will  also  remain  the  same  as  in  Example  1.  The  displace¬ 
ment  response  for  the  model  and  the  structure  can  be  seen 
in  Figures  20-29.  As  in  the  last  exao^le  both  the  model 
response  and  the  struct\xre  response,  that  includes  the 
residuals,  move  towards  the  Ideal  model  as  the  number  of 
Inputs  increase.  At  K:9  the  structure  response  is  again 
nearly  identical  to  the  response  of  the  ideal  model. 


Table  2.  Example  2  closed-loop  eigenvalues 


IDEAL 

i::2 

i::9 

-1.  23 

-1.  44 

-1.  23 

-78.  76 

-75.  29 

-78.  76 

-33.  56 

-3.  43-137.  79 

-33.  64 

-46.  43 

-3.  43+137.  79 

-46.  32 

-40.  00*179.  31 

-6.  761188.  57 

-39.  92*179.  35 

-40.  00*1152.  76 

-9.  11*1157.  14 

-39.  66*1152.  85 

-40.  00*1243.  48 

-8.  89*1245.  23 

-39.  11*1243.  62 

-40.  00ii353.  05 


-6.  63*1353.  63 


-38.  75*1353.  19 


If  we  examine  Table  2  we  can  compare  tbe  closed- loop 
eigenvalues  of  tbe  Ideal  model  to  tbat  of  tbe  actual  model 
at  K=2  and  at  Ts9.  From  tbls  table  we  can  confirm  tbat  tbe 
actual  model  is  almost  at  tbe  ideal  case  wben  K=9.  Tbe 
input  force  response  for  K=2-9  can  be  seen  in  Figures  30-33 
and  we  can  see  tbat  as  tbe  number  of  inputs  increase  tbe 
input  force  decreases  towards  tbe  ideal  distributed  case. 


Pig.  21.  Example  2  structure  displacement/time  K=2 


Fig.  23.  Example  2  structure  displacement /time  K=4 


Fig.  24.  Example  2  model  displacement/time  K:6 


Fig.  25.  Example  2  structure  displacement/time  K:6 


In  ttxls  example  the  control  system  will  he  applied  to 
basically  the  same  simply  supported  beam,  but  the  length 
will  be  increased  to  L=10  units.  All  the  other  physical 
qualities  of  the  beam  will  remain  the  same  as  in  the 
earlier  examples.  The  structure  will  consist  of  Ns=10 
modes  and  the  model  will  have  N=6  modes.  The  damping 
ratio  is  a:~4.  0  and  the  stiffness  coefficient  p=0.  The 
observed  point  is  selected  as  x=2.  1  units.  The  time  tics 
on  the  plots  are  at  4  second  intervals  and  the  displacement 
and  force  tics  are  at  0.2  unit  intervals.  When  Zz6  the 
input  points  are  at  x:0.  9,  2.1,  3.1,  4.2,  5.3,  6.1,  7.3, 
and  8.2  units.  When  K=10  the  points  x=9.  3,  9.7  are  added 
and  for  K:14  the  points  x=4.  75,  1.6,  0.5,  and  3.55  units 
are  included.  When  K:17  the  final  input  points  are 
x=6. 45,  2.  7  and  O.  7  units.  The  displacement  responses  for 
the  model  and  the  structure  are  shown  in  Figures  34-41. 

By  exauninlng  the  responses  we  note  that  by  K=8  the  model 
appears  to  have  already  reached  the  ideal  case.  If  we 
examining  Table  3  we  can  see  that  the  closed- loop  eigen¬ 
values  are  not  close  enough  to  the  ideal  eigenvalues  and 
the  structure  still  has  a  very  large  residual  effect. 
However,  by  K=17  the  structure  response  closely  resembles 
the  ideal  model  with  only  a  small  amount  of  residual 


remaining. 


Table  3.  Example  3  closed-loop  eigenvalues 


IDEAL 

K=8 

K:17 

0024367 

-.  0024368 

-. 0024367 

039351 

-.  039388 

-. 039356 

20806 

20968 

-.  20827 

77267 

81122 

-. 776622 

-3.  2273 

-3.  0980 

-3.  2110 

-3.  7919 

-3.  7699 

-3.  7885 

-3.  9606 

-3.  9577 

-3.  9601 

-3.  9976 

-3.  9974 

-3.  9975 

-2.  00111.  445 

-1.  85111.  659 

-1.  98111.  464 

-2. 00112.  937 

-1.  61113.  119 

-1.  97112.  956 

CHAPTER  in 


CONTROL  OP  AN  UNDAMPED  TRUSS 


Introduction 

In  this  chapter  the  closed- loop  feedback  control 
method  and  the  open- loop  method  will  be  applied  to  a  finite 
degrees  of  freedom  undanqped  truss  structure.  We  will 
examine  the  structure  and  the  implementation  of  the  control 
system  on  this  structure,  then  we  will  look  at  two  closed- 
loop  and  two  open- loop  examples. 

Implementation 

The  truss  that  we  will  be  applying  the  control  system 
to  [l£]  is  a  two-dimensional  undamped  structure  (Figure 
42)  composed  of  seven  beams  with  unit  mass  distribution, 
unit  stiffness  EIA=l,  and  unit  length  L=l.  The  struct\ires 
left  side  is  pinned  to  a  wall  at  both  Joints.  The  lower 
left  Joint  where  beams  1  and  2  are  connected  is  the 
origin  for  the  global  coordinates  x  and  y.  Ihe  generalized 
coordinates  used  to  describe  the  displacement  at  each 
Joint  are  Uj,  U2,...Uio»  but  since  two  of  the  joints  are 
pinned  their  displacements  will  be  zero  and  their  general¬ 
ized  coordinates  are  dropped  from  the  final  equation. 
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(76 


nie  equation  of  motion  of  the  truss  has  the  form 
mi  +  KU  =  DF 

where  H  Is  the  mass  matrlXi  K  Is  the  stiffness  matrix,  F 
Is  the  Input  force  vector,  and  D  Is  a  force  distribution 
matrix.  The  mass  and  stiffness  matrix  have  the  form 
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-1 

0  0 

-1 

1 

0 

l+2f2 
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1 

-1 

0 

0 
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-2f2 
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The  displacement  and 
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force  vector 
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form 
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The  closed-loop  control  on  the  truss  Is  liqpleroented 


In  the  same  manner  as  on  the  beam  with  the  exception  that 
the  Input  forces  for  these  example  are  only  being  applied 
at  the  Joints  and  In  the  direction  of  the  structure's 
degrees  of  freedom  as  depicted  in  Figure  42.  The  computer 
simulation  of  the  closed-loop  system  Is  accomplished 
through  the  program  PHTCHTF  (Appendix  B). 

The  open-loop  control  is  accomplished  by  controlling 
the  model  without  any  feedback  from  the  structure  on  the 
effect  of  the  residual  modes.  When  the  open- loop  control 
method  Is  Implemented  the  model  response  for  any  number  of 


input  points  will  always  be  the  ideal  case,  but  the 
structure  response  will  move  towards  the  ideal  model  as 
the  number  of  inputs  increase.  The  conqputer  simulation  of 
the  open- loop  system  is  accomplished  by  modifying  the 
program  PHTCNTF  so  that  there  is  no  feedback  from  the 
structure  to  the  model. 

For  all  the  examples  the  initial  modal  state  distur¬ 
bances  are  chosen  as  qp(0)=0  and  qp(0)=1.0  units/sec, 
r=l,2,  ...N.  The  Joint  chosen  to  observe  the  model  and 
structure  responses  will  be  the  upper  right  Joint  where 
the  beams  4,5,  and  7  are  connected.  Since  the 
structure  is  two  dimensional  we  will  examine  the  displace¬ 
ment  and  input  force  in  both  the  x  and  the  y  direction. 
The  horizontal  axis  on  all  the  plots  will  reflect  the  time 
and  the  tics  at  2  second  intervals.  The  vertical  axis 
will  reflect  either  the  displacement  or  the  input  force 
with  the  tics  at  0.25  unit  intervals.  The  structure  has  a 
finite  number  of  modes,  Ns=6,  and  the  reduced  order  model 
controls  N:n  modes  where  N<Ng.  The  number  of  inputs  will 
be  referred  to  as  E. 

Example  4 

In  this  closed-loop  control  system  example,  the 
reduced  order  model  will  consist  of  the  first  two  modes  of 


the  structure  N=2.  The  damping  coefficient  is  selected  as 
a=-2  and  the  stiffness  coefficient  is  selected  as  E=0. 


When  the  Input  forces  P5  and  P5  are  being  applied. 

When  the  input  forces  P7  and  Pg  are  added  and  for  K=6 

the  input  forces  P9  and  P^q  are  included.  The  model  and 

structure  displacement  responses  are  shown  in  Pigures 

43-54.  As  in  the  previous  examples  as  the  number  of 

inputs  increase  both  the  actual  model  and  Y.he  structure 

approach  the  ideal  model  response  which  is  reached  when 

E=6.  If  we  look  at  the  ideal  gain  matrix  for  control  in 

the  X  direction  we  see  that  it  has  the  form 

'0  0  0  o' 

0-200 

G  =  0  O  0  O  (81) 

0  0  0  -2  J. 

When  k;=2  the  non-ideal  gain  matrix  is 

■  0  0  0  o' 

0  -0.  07  0  -0.  05 

G  =  0  0  0  0  (82) 

”  0  -O.  05  0  -0.  04  . 

When  K=4  the  non-ideal  gain  matrix  is 

0  0  0  o' 

0  -1.20  0  -0.59 

G  =  0  0  0  0  (83) 

^  0  -0.59  O  -O.  69J. 

When  K:6  the  non-ideal  gain  matrix  is 

■  0  0  0  o' 

O  -2.  0  0  6.  IxlO-l'^ 

G  =  0  0  0  0  (84) 

”  [  O  6.  IxlO-i^  0  -2.0 

By  examining  the  gain  matrices  it  can  be  easily  seen  that 

as  the  number  of  inputs  increase  the  actual  model,  that 

includes  feedback  from  the  structures  residual  modes, 


approaches  the  ideal  model.  By  examining  Table  4  we  can 


nils  closed-loop  exaoqple  Is  very  similar  to  Ezanqple  4 
except  tnat  the  reduced  order  model  Is  Increased  to  Include 
the  first  three  modes  of  the  structure  H=3  and  the  damping 
coefficient  is  Increased  to  0(=-5.  All  the  other  param¬ 
eters  remain  the  same  as  in  Example  4  including  the  input 
points.  The  displacement  response  for  the  model  and  the 
stznicture  in  both  the  x  and  y  direction  can  be  seen  in 


Figures  57-68. 


In  ttiis  example  we  will  apply  ttie  open- loop  control 
system  to  tlie  undamped  truss.  To  get  a  good  comparison  of 
the  open-loop  method  and  closed-loop  method  we  will  apply 
the  open- loop  control  to  the  same  problem  as  In  Exanqple  4. 
Therefore,  the  reduced  order  model  will  be  H=2  and  the 
damping  and  stiffness  coefficients  will  be  a=-2  and  3=0. 

Ihe  Inputs  will  be  the  same  as  In  Example  4.  When  open- 
loop  control  Is  used  the  effect  of  the  residual  modes  on 
the  structure  Is  not  feed  back  to  the  reference  controlled 
model,  equation  (32).  Therefore,  the  reference  model 
response  will  always  reflect  the  Ideal  case,  this  can  be 
seen  In  Figures  69  and  73.  The  structure  displacement 
responses  In  both  the  x  and  y  direction  can  be  seen  In 
Figures  70-72  and  74-76.  In  examining  these  figures  we  can 
see  that  the  effect  of  the  open- loop  control  on  the 
structure  Is  not  as  strong  as  the  closedloop  control  when 
the  Input  number  Is  small.  However,  when  enough  Input 
points  are  used  the  structure  will  resemble  the  Ideal  model 
and  the  effect  of  the  residual  modes  is  eliminated.  In 
Figures  77  and  76  we  can  see  that  the  Input  force  response 
decreases  as  the  n\imber  of  Inputs  Increase. 


73 

Example  7 

In  this  example  we  will  apply  the  open- loop  control 
system  to  the  same  problem  as  in  Example  5.  The  reduced 
order  model  will  consist  of  the  first  three  modes  of  the 
structure  N=3  and  the  damping  and  stiffness  coefficients 
will  be  a=-5  and  p:0.  All  the  input  points  remain  the 
same  as  in  Example  5.  The  ideal  model  displacement 
response  in  both  the  x  and  y  direction  can  be  seen  in 
Figures  79  and  83.  The  structure  displacement  responses  in 
both  the  X  and  y  direction  can  be  seen  in  Figures  80-82 


and  84-86. 


CHAPTER  IV 


COHCLUSIOHS 

In  this  thesis  we  have  heen  examining  a  new  approach 
forwarded  in  [2]  for  controlling  large  flexible  structures. 
The  importance  of  this  new  approach  is  that  by  using  the 
independent  modal -space  control  method  and  a  finite  number 
of  spatially  distributed  input  points  one  can  essentially 
eliminate  control  spillover  into  modes  beyond  the  reduced 
order  model  of  the  system.  By  examining  the  formulation  of 
this  control  method  we  see  that  the  control  of  the  system 
is  accomplished  through  designing  the  input  control  forces 
in  the  form  of  equation  (56).  This  spatial  modal  Input- 
distribution  control  redistributes  the  spatially  discrete 
inputs  so  as  to  create  effective  continuously  distributed 
damping  and  stiffness  inputs.  The  inherent  spatial 
orthogonal  non-pass  filtering  property  of  this  control 
results  in  the  essential  elimination  of  the  excitation 
effect  of  the  residual  modes. 

In  Cl:iapter  II  and  III  we  applied  the  closed- loop 
control  method  to  an  undamped  simply  supported  beam  and  an 
undamped  timiss.  By  examining  the  displacement  responses 
of  the  reduced  order  model  and  the  structure  we  can  see 


that  as  the  number  of  Inputs  increase  the  reduced  order 
model  approached  the  Ideal  model  case  where  the  effect  of 
residual  mode  excitation  is  eliminated  from  the  structure. 
If  we  compare  the  reduced  order  model's  gain  matrix  and 
closed-loop  eigenvalues  to  those  of  the  ideal  model  we  can 
see  that  the  reduced  order  model  approaches  the  ideal 
model  with  a  finite  number  of  inputs.  By  looking  at  the 
input  force  responses  we  can  see  that  the  input  force  is 
reduced  with  the  increase  of  inputs  until  it  reaches  the 
ideal  input  force  response  which  represents  a  continuously 
distributed  input  force.  We  also  applied  an  open- loop 
control  system  to  the  truss  problem  and  by  examining  the 
displacement  responses  we  can  see  that  with  enough  input 
points  the  residual  mode  excitation  can  be  eliminated  Just 
as  effectively  as  the  closed-loop  control  system.  It  is 
very  important  to  note  that  higher  numbers  of  input  points 
were  examined  for  each  of  the  examples  shown,  but  the 
effect  on  the  response  plots  is  so  small  that  it  is 
Impossible  to  see. 

Prom  these  examples  it  is  easy  to  see  that  by  using 
spatial  modal  input-distribution  control  and  a  finite 
number  of  inputs  it  is  possible  to  essentially  eliminate 
control  spillover.  Therefore,  the  method  forwarded  in  [2] 
can  be  a  very  useful  approach  in  controlling  the  large 
flexible  space  structures  of  the  future. 


//  JOB  , 

//  TIHE= (2,0),REG10N=2560t 
/«JOBPARM  LINES=9000 
//  EXEC  PLOTV77 

//SYSLIB  DD  DSNrSYSl.  VFORTLIB,  DISPrSHR 
//  DD  DSNrSYSl.  PLOTLIB,  DISP=SHR 

//  DD  DSH=SYS1.  PORTLIB,  DISP  =  SHR 

//  DD  DSN=SYSl.  IHSL.  DOUBLE,  DISP=SHR 

//GO.  SOURCE  DD  « 


PROGRAM  PHTCNT 


IMPLICIT  REAL»ia(A-H,  0-Z) 

EXTERNAL  F1.F2 
DOUBLE  PRECISION  CF 

DIMENSION  Cl  1 (40,  40) ,  C (40,  40) ,  CCDUM(40,  40) ,  WCS (40) , 
«CCINV(40,  40) ,  WKAREA(500) ,  A(40,  40) ,  BCB2 (40,  40) , 
«WC(40) ,  WSTS (40) ,  WSTINS (40) ,  GM(40,  40) ,  PHIM(40,  40) , 
»G(40,  40) ,  ALPHA (40),  OMEGA (40) ,  BETA (40) ,  B(40,  40) , 

» GAMS (40,  40) ,  GAMH (40,  40) ,  BS (40,  40) ,  BCB (40,  40) , CBW(40) 
«PHIS(40,  40) ,  BT(40,  40) ,  CINVBT(40,  40) ,  WSTINH(40) , 
«GHN(40,  40) ,  ACLI (40, 40) ,  ACLN(40, 40) , WK (1700) , WSTH(40) 
COMPLEX# 16  EVALI (40) ,  EVALN (40) , EYECl (40, 40) , 
mZNI, ZNN, EVECM(40,  40) 

COMMON  X (40) ,  STIFF,  DLEN,  DNORM,  N,  KPM,  H,  IDUH,  DUHJ,  PI 
COMHON/C0MCIZ/DISP(5002,  46) ,  YM ( 10) ,  YT ( 10) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c  c 

C  INPUT  SYSTEM  INITIAL  CONDITIONS  C 
C  C 

ccccccccccccccccccccccccccccccccccccccc 

READ#,  N,  NS,  K 

READ# ,  DT,  STIFF,  DLEN,  RO,  NUM 
READ#,  ALPHA1,BETA1 
DO  10  1  =  1, K 
10  READ#,  X(I) 

N2:N#2 
NS2=NS#2 
LD=40 
NDIM=2#K;4  1 

DN0RM=DSQRT(2.  DO/DLEK) 

PI  =  DATAN(1.D  00)  #4. 

IDGT=8 

IJ0B=2 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  INITIAL  CONDITIONS  FOR  PLOT  C 
C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
YH(1)  =  30. 


YH(2)  =  30. 

YH(3)  =  30. 

YM(4)  =30. 

YT(1)  =  10. 

YT(2)  =  10. 

YT(3)  =  10. 

YT(4)  =  10. 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  INPUT  EIGENVALUES  C 

C  C 

ccccccccccccccccccccccccccccccccc 

DO  15  1=1. NS 
RI  =  I 

15  OHEGA(I)  =RlKK2KPlMK2iiDSQRT(STIPF/RO)/DLENH«2 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C  C 

C  IHSC  MODAL  GAIN  MATRIX  C 

C  C 

ccccccccccccccccccccccccccccccccccc 

DO  16  1  =  1. N 

16  BETA(I):BETA1/0MEGA(I) 

DO  17  1=1, N 

17  ALPHA(I) =ALPHA1 
DO  20  1:1.  N2 
DO  20  J:1,N2 

20  QM(I,  J)=0, 

DO  25  1  =  1,  N 
GM(2»1,  2«I-1)  sBETAd) 

25  GM(2»I,  2x1) :ALPHA(I) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  INPUT  MODAL  MODEL  INITIAL  CONDITIONS  C 
C  C 

cccccccccccccccccccccccccccccccccccccccccc 

DO  151  1  =  1, H2 
READ*,  WSTINH(I) 

151  WSTM(I) =WSTINH(I) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c  c 

C  STRUCTURE  INITIAL  CONDITIONS  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  145  1=1, NS2 
WSTS(I)  =0. 

145  WSTINS(I)=0. 

DO  146  1=1, N2 
WSTINS(I) =WSTINM(I) 

146  WSTS(I) =WSTINS(I) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  C  MATRIX  C 
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cccccccccccccccccccccccccccc 
DO  4  1:1,  K 
DO  4  J:1,K 

IF(X(I).  LE.  X(J)  )  C(I,  J)  :CF(X(I),  X(J)  ) 

4  IF(X(I).  GT.  X(J)  )  C(I,  J)=CF(X(J),  X(I)  ) 

DO  49  I:1.E 
DO  49  J:i,E 
49  011  (1,  J):C(1,  J) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  B.  BT,  AND  BS  MATRICES  C 
C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  12  1:1. N 
DO  12  J:1,K 
DI:I 

B(I,  J) =DNORMhDSIH(DImX(J) hPI/DLEK) /0MEGA( I ) 

12  BT(J,  I):B(I,  J) 

DO  13  I:1.HS 
DO  13  J:1,K 
DI:I 

13  BS(I.  J)  :DHORHnDSIN(DIiiX(J)  mPI/DLEN) /OMEGA ( I ) 

ccccccccccccccccccccccccccccccccccc 
c  c 

c  C  IHV.  AND  BCB  MATRICES  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CALL  LIHV2F(C,  K,  LD.CCINV,  IDGT.WKAREA,  lER) 

CALL  GPROD(CCIHV,  BT,  CINVBT,  K,  K,  N) 

CALL  GPROD (BS, CINVBT,  BCB,  NS,  K,  N) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  TRANSITION  &  GAMMA  MATRICES  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC  ccccc 
DO  500  1:1, NS2 
DO  500  J:1,NS2 
GAMS  (I,  J)  :0. 

500  PHIS(I,J):0. 

DO  501  1:1,  N2 
DO  501  J:1,N2 
GAMHd,  J):0. 

501  PHIM(I,J):0. 

DO  502  1:1, NS 

PHIS(2iiI-l,  2«I-1)  :DCOS(OMEGA(I)  »DT) 

PHIS(2nI-1,  SkI) :DSIN(0MEGA(I) kDT) 

PHIS(2»I,  2»I-1) :-PHlS(2»I-l,  2«I) 

PHIS (2k 1,  2kI) =PH1S(2kI-1,  2kI-1) 

GAMS(2kI-1, 2HI-1) :DSIH(OMEGA(I) kDT)/OMEGA(I) 
GAMS(2kI-1,  2kI) : (1.  -DCOS(OHEGA(I) kDT) )/OHEGA(I) 
GAMS  (2k  I, 2kI-1)=-GAMS(2kI-1,  2k  I) 


502  GAHS(2iiI,  2)il)  =aAHS(2iiI-l.  2iil-i) 

DO  503  1=1, N2 

DO  503  J=1,N2 
PHIHd,  J)=PHIS(1,  J) 

503  GAHHd,  J)  =aAHS(I,  J) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c  c 

C  CLOSED  LOOP  EIGEHVALUES  &  VECTORS  C 
C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  201  1:1,  H2 
DO  201  J=1,H2 
Ad,  J)  =0. 

ACLI  (I,  J)=0. 

ACLMd,  J)  =0. 

GHNd,  J)  =.  0 

201  BCB2d,J)=0. 

DO  202  1=1, N 
DO  202  J:1,N 

202  BCB2(2mI,  2iiJ)  =BCBd,  J) 

CALL  GPR0D(BCB2,GH,  GHN,  N2.H2,H2) 

DO  204  1:1, N 
A(2iil-1,  2iil)  :OHEGAd) 

204  A(2iil,  2iil-l)  :-OHEGAd) 

CALL  GHADD ( A,  GH,  ACL I ,  N2 ,  N2 ) 

CALL  GHADD (A,  GHH,  ACLH,  K2,  H2 ) 

CALL  EIGRP(ACLI,  H2,  LD,  IJOB,  EVALI,  EVECI,  LD.  WE,  lER) 
CALL  EIGRF(ACLN.  N2,  LD,  IJOB,  BVALN,  EVECN,  LD,  WE,  lER) 
CCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  PRINT  OUTPUT  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCC 
VRITE(6,  301)  N,  NS,  E 

301  FORSAT(//,  2X,  'N=',  14,  4X,  'NS=',  14,  4X,  'E=',  14) 
WRITE(6,  302)  DT,  STIFF,  DLEN,  RO,  NUM 

302  FORWAT{//,  2X,  'DT=*,  E12.  4,  2X,  'STIFF:  ' ,  E12.  4,  2X, 
n'DLEN: ',E12.  4,  2X,  'R0:',E12.  4,  EX,  'NUH=',  13) 

WRITE (6, 303)  ALPHAl,  BETAl 

303  FORMAT  (//,  2X,  'ALPHA:  ' ,  E12.  4,  4X,  '  BETA:  ',E12.4) 
WRITE (6, 304) 

304  FORMAT (//,  2X,  '  X  VECTOR  ') 

WRITE(6,  100)  (Xd),  1:1, E) 

WRITE(6, 305) 

305  FORMAT (//,  2X,  '  WSTINM  VECTOR  ') 

WRITE(6,  100)  (WSTINHd),  1=1,  N2) 

100  FORHAT(//,  2X,  7(2X,  D16.  7) ) 

WRITE(6,  107) 

107  F0RHAT(//2X,  'C  MATRIX') 

DO  21  1=1,  E 

21  WRITE(6,  100)  (Cll (I,  J),  J=l,  E) 

WRITE(6, 111) 


o  o  o  o 
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111  FORMAT (//,  2X,  'C  INVERSE  MATRIX' ) 

DO  40  I:1,K 

40  WRITE (6,  100)  (CCINV(I,  J) ,  J: 1,  K) 

WRITE(6.  152) 

152  FORMAT (//,  2X,  '  BOB  MATRIX') 

DO  153  1=1. NS 

153  WRITE(6.  100)  (BOB ( I,  J) ,  J= 1 ,  N) 

WRITE(6,  170) 

170  FORMAT{//,  2X,  'IDEAL  GAIN  MATRIX') 

DO  171  1  =  1, H2 

171  WRITE(6.  100)  (GM{I,  J),  J=1,N2) 

WRITE(6,  172) 

172  FORMAT (//,  2X,  'NON- IDEAL  GAIN  MATRIX') 

DO  173  1=1, N2 

173  WR1TE(6,  100)  (GMN(I,  J) ,  J= 1,  N2) 

WR1TE(6, 161) 

161  FORMAT (//,  2X,  'COMPLEX  CLOSED  LOOP  EIGENVALUES') 
WRITE(6,  160)  (EVALN(I),  1=1,  N2) 

160  FORMAT(//,  2X,  3(2X,  2E16.  7)  ) 

WRITE(6, 162) 

162  FORMAT(//,  2X,  'IDEAL  COMPLEX  CLOSED  LOOP 
K  EIGENVALUES') 

WRITE(6,  160)  (EVALI (I),  1:1,  N2) 

WRITE{6,  163) 

163  FORMAT (//,  2X,  'COMPLEX  CLOSED  LOOP  MODAL  MATRIX') 
DO  164  1=1, N2 

164  WRITE (6,  160)  (EVECN(I,  J) ,  J= 1,  N2) 

WRITE(6,  165) 

165  FORMAT (//,  2X,  'IDEAL  COMPLEX  CLOSED  LOOP  MODAL 
«  MATRIX') 

DO  166  1  =  1,  N2 

166  WRITE(6,  160)  (EVECI (I,  J) ,  J= 1,  N2) 

WRITE(6,  150) 

150  PORMAT(//,  2X,  'PHI  MATRIX'  ) 

DO  154  J=1,NS2 

154  WRITE{6,  100)  (PHIS(I,  J),  1=1,  NS2) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

STRUCTURAL  RESPONSE  TO  MODAL  INPUTS  C 

(WSTS= STRUCTURE  MODAL  COORDINATES)  C 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
TIME: -DT 
KTIMErO 

147  KTIHE:KTIME+1 
TIME:TIME+DT 

IF  (KTIME.  EQ.  1)  GO  TO  510 
DO  136  1=1, N2 
WC(I)  =0. 

DO  138  J=1,N2 

138  WC(I)  =WC{I)  ♦GMd,  J)  «WSTS(J) 

DO  504  1=1, NS2 
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504  WCS(I) sO. 

DO  505  1:1. NS 
WCS(SkI)=0. 

DO  505  J:1,H 
CBW(J)  :0. 

505  VCS(SiiI)  :VCS(2«1)+BCB(I.  J)iiWC(2iiJ) 

DO  141  1:1. NS2 

WSTS(I)  :0. 

DO  141  J:1.HS2 

141  VSTS(I):WSTS(I)-^PHIS(I.  J)kVSTIHS(J) 

DO  142  I:1.NS2 

DO  142  J:1.NS2 

142  WSTS  ( I )  :  VSTS  ( I )  4<3AHS  ( I.  J)  kWCS  (  J) 

DO  143  1=1. HS2 

143  WSTIHS(I) =VSTS(1) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c  c 

C  MODEL  RESPONSE  TO  MODAL  INPUTS  C 

C  (VSTM= MODEL  MODAL  COORDINATES)  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  136  I:1.N2 
WSTM(I)  =0. 

DO  136  J:1.N2 

136  WSTM(I)  =WSTM(I)+PHIM(I.  J)«WSTINM(J) 

DO  506  1:1.  N2 

DO  506  J=1,N2 

506  WSTM(I)sWSTM(I)+aAMM(I.  J)iiWCS(J) 

DO  137  1:1.  N2 

137  WSTINM(I) :WSTM(I) 

DO  160  1=1, K 

DO  160  J=1,N 

160  CBW(I)  :CBW(I)  +CINVBT(I.  J)  )iWC(2iiJ) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  MODEL  &  STRUCTURE  C 

C  DISPLACEMENT  MATRIX  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
510  DO  506  J=1.E 

DISP(KTIME.  J*l)  =  0. 

DISP(KTIME.  J+l+K)  =0. 

DISP(KTIME,  J+l+i:+K)=0. 

DO  507  1:1,  N 
DI  =  I 

507  DISP(KTIME.  1 )  :DISP  (CTIME.  1 )  ^^DNORMnDSIN  (DI  hX  ( J) 

imPI/DLEN)  iiWSTM(2mI-1) 

DISP(KTIME.  J+l+K+K)  =CBW(J) 

DO  509  1  =  1.  NS 
D1  =  I 

509  DISP(KTIME,  J+l+K)  =DISP(KTIME,  J+l+K)  +DNORMi»DSIN 
N  (DInX(J)  nPI/DLEN)  iiWSTS(2iiI-1) 


DISP  (KTIHE,  Jf  1  -fE)  =  10.  aDISP  (mHE,  1 ) 

DISP (ETIHE,  J+ 1 )  =  10.  aDISP (ETIHE,  J+ 1 ) 

DISP (ETIHE,  J+ 1 ♦E+E)  =  1 0.  aDISP (ETIHE,  J+ 1 ♦E+E) 

508  CONTINUE 

DISP (ETIHE,  1)=TIHE 

IF(TIHE.  LT.  2.  5D  00)  GO  TO  147 

WRITE(6,  70) 

70  FORHAT(//,  2X,  'CIHVBT  HATRIX' ) 

DO  71  1  =  1,  E 

71  WRITE(6,  100)  (CINVBTd,  J),  J=l,  N) 

CALL  HPLOT(2.  5D  00,0.  5D  00,  ETIHE,  E,  HUH) 

STOP 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  INFLUENCE  FUNCTION  (CF)  C 

C  C 

ccccccccccccccccccccccccccccccc 

DOUBLE  PRECISION  FUNCTION  CF(Y,  ETA) 

IHPLICIT  REALaa(A-H,  0-2) 

COHHON  X (40) ,  STIFF,  DLEN,  DNORH,  N,  E,  H,  IDUH,  DUHJ,  PI 
CF:  (Yaa3aETA/DLEN>YaETAaa3/DLEH-f2.  D  OOaYaETAaDLEN 
a-Yaa3-3.  D  00aYaETAaa2)/6.  D  OO/STIFF 
RETURN 
END 

CCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  HATRIX  PRODUCT  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  GPROD (A,  B,  C,  H,  N,  E) 

IHPLICIT  REALa8(A-H,  0-2) 

DIHENSION  A (40,  40) ,  B(40, 40) , C(40, 40) 

DO  1  1  =  1, H 
DO  1  J:1,E 
C(I,  J)  =0. 

DO  1  L=1,N 

1  C(I,  J)  =C(I,  J) +A(I,  L)  aB(L,  J) 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  HATRIX  ADDITION  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  GKADD(A,  B,  C,  N,  H) 

IHPLICIT  REALa6(A-H,  O-Z) 

DIHENSION  A (40,  40) ,  B(40,  40) ,  C (40,  40) 

DO  1  1  =  1, N 
DO  1  J=1,H 

1  C(I,  J)  =A(I,  J) +B(I,  J) 

RETURN 


ccccccccccccccccccccccccc 
c  c 

C  PLOTTING  ROUTINE  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  HPLOT (XHAG. XT,  H,  KH. KPT) 

IMPLICIT  REALii8(A-H.  0-Z) 

COMHON/COHC I Z/A ( 500S ,  46 ) ,  YM ( 1 0 ) , YT ( 1 0 ) 

XSCALE=10.  /2.  5 

XORQ=l. 

YORG=5.  5 

CALL  PLOTS  (0,0,0) 

DO  5  1=1,3 

CALL  PLOT(l.  0,  0.  ,  -3) 

CALL  PLOT(0.  ,  11.  ,  2) 

CALL  AX(5.  OD  00,  1.  OD  00,  YH  ( I ) ,  YT  ( I ) ,  XORG,  YORG) 
YSCALE=3.  5/YH(I) 

K=3 

IPd.EQ.  1)  NrKPT^l 
IFd.EQ.  2)  N=KH+KPT+1 
IFd.EQ.  3)  N:KH+EH4^KPT^1 
DO  1  J=1,M 
X  =  A(J,  1)«XSCALE+X0RG 
Y=A(J,  N) nYSCALE+YORG 
CALL  PLOT(X,Y,K) 

1  CONTINUE 

CALL  PLOTdO.  5,  0.  , -3) 

5  CONTINUE 

CALL  PL0T{0.  ,  0.  ,  999) 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  PLOTTING  SUBROUTINE  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  AX (XHAG,  XT,  YMAG,  YT,  XORG,  YORG) 
IMPLICIT  REAL«8(A-H,  O-Z) 

XH=10. 

YM:3.  5 

CALL  PLOT(XORG,  YORG-YM-.  5,  3) 

CALL  PLOT (XORG,  YORG +YH+.  5,  2) 

Y:YM«YT/YHAG 
YY  =  YH 
YI  =  YH+.  1 

1  CALL  PLOT(XORG,  YORG+YY,  3) 

CALL  PLOT (XORG-.  1, YORG+YY,  2) 

YY=YY-Y 

IF(YY.  GE.  -YI)  GO  TO  1 
X  =  XHiiXT/XHAG 
XX  =  XH 


CALL  PLOT ( XORG,  YORG,  3 ) 

CALL  PLOT ( XORG ♦XH+.  5.  YORG, 2) 
2  CALL  PLOT(XORG+XX,  YORG,  3) 
CALL  PLOT(XORG+XX,  YORG-.  I, 2) 

xx=xx-x 

IF(XX.  GT.  .  1)  GO  TO  2 

RETURN 

END 

/» 

//GO.  DATA  DD  « 

(  Place  Data  Here  ) 


/» 

//GO.  PLOTPARH  DD  « 
&PLOT  XHAX=120  &END 


//  JOB  , 

//  TIHE= (2.  0).  REGI0H=2816C 
/kJOBPARH  LINES: 9000 
//  EXEC  PLOTV77 

//SYSLIB  DD  DSH=SYS1.  VFORTLIB,  DISP=SHR 
//  DD  DSNrSYSl.  PLOTLIB.  DISP=SHR 

//  DD  DSHsSYSl.  FORTLIB,  DISPrSHR 

//  DD  DSN=SYS1.  IHSL.  DOUBLE.  DISP:SHR 

//GO.  SOURCE  DD  m 


PROGRAM  PHTCNTF 
(FINITE  MODEL) 


IMPLICIT  REALk8(A-H.O-Z) 

DIMENSION  Cll (40,  40).  BSS(40,  40),  EV(40.  40),WCS(40}, 
«CCIHV(40,  40} ,  WKAREA (500) ,  A (40,  40) .  BCB2 (40,  40) , 
kWC  (40) . VSTS (40) , VSTINS(40) . GN(40.  40) .  PHIH(40,  40) , 
me (40,  40) .  ALPHA (40)  ,  OMEGA (40) , BETA (40) ,  B (40. 40) . 
mGAMS(40,  40) ,  GAMM(40,  40)  ,  BS(40.  40) ,  BCB(40.  40) , 
mPHIS (40,  40) ,  BT(40,  40) , CINVBT(40,  40) ,  WSTINM (40) , 
mGMH (40.  40) .  ACLI (40,  40) ,  ACLN (40,  40) ,  WE (1700) . 
mSTIFM (40,  40) ,  DMASS (40, 40) ,  DK (40,  40) ,  DH(40,  40) , 
hDDK(40,  40) ,  DDH(40,  40) ,  D(40,  40) ,  DHI (40,  40) ,  CBW(40) , 
mAOL(40, 40) .  0MEGSQ(40) , BT(40.  40) .  DMEV(40,  40) . 
hWSTH(40)  .  DHM(40) ,  DMORH(40) 

COHPLEXm 16  EVALl (40) ,  EVALN (40) ,  EVECI (40,  40) , 
hZNI,  ZNN,  EVECM(40,  40)  .  CW(40)  .  CEV(40,  40) 
COMMOB/COHCIZ/DISP(5O02,  46) ,  YM(IO) ,  YT(IO) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c  c 

C  INPUT  SYSTEM  INITIAL  CONDITIONS  C 
C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
READm,  N.  NS.K 

READm,  DT,  STIFF,  DLEN,  RO,  NUH 
READm,  ALPHAl.BETAl 
DO  7  1=1, NS 

7  READm,  (D(I,  J) ,  J:1.K) 

DO  5  1:1,  NS 

5  READm,  (STIFHd,  J)  ,  J:  1,  NS) 

DO  6  1=1, NS 

6  READm,  (DMASS ( I.  J) .  J: 1 .  NS) 

N2=Nh2 

NS2=NSm2 

LD=40 

KDIM=2mK41 

IJ0B=2 

IDGT:6 

IJOB:2 

CZ:STIFF/(2mDSQRT(2.  D  OOmDLEN) ) 


CH=R0nDLENmii2/(6mDSQRT(2.  D  00)) 

DO  9  1=1. HS 
DO  9  J=1.NS 
DK(I,  J)=CZi»STIPH(I,  J) 

DH(I,  J)=CHmDHASS(1.  J) 

DDK(I.  J)=DK(I.  J) 

9  DDH(I.  J)  =DH(I.  J) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  INITIAL  CONDITIONS  FOR  PLOT  C 
C  C 

ccccccccccccccccccccccccccccccccc 

yM(l)  =  30. 

YH{2)  =  30. 

YH(3)  =  30. 

YT(1)  =  10. 

YT(2)  =  10. 

YT(3)=10. 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  INPUT  MODAL  MODEL  INITIAL  CONDITIONS  C 
C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  12  1=1. H2 
READn.  WSTIHH(I) 

12  WSTM(I)  =VSTIHN(I) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c  c 

C  STRUCTURE  INITIAL  CONDITIONS  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  13  1=1. HS2 
WSTS(I)  =0, 

13  WSTINS(I)=0. 

DO  14  1=1, N2 
WSTINS ( I ) =WSTINH ( I ) 

14  WSTS ( I ) =WSTINS (I ) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c  c 

C  FIND  EIGENVALUES  &  VECTORS  C 
C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CALL  LINV2F(DN.  NS,  LD,  DMI,  IDGT.  WKAREA,  lER) 
CALL  GPROD(DMI,DK,  AOL.  NS,  NS.  NS) 

CALL  EIGRF(AOL.  NS,  LD,  IJOB,  CV,  CEV,  LD,  WE,  lER) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  NORMALIZED  MODAL  MATRIX  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  20  1=1, NS 
DNORM(I)  =  0. 


OHEGSQ ( I ) =RBAL (CW ( I ) ) 

OMEGA ( I ) =  DSQRT (OHEGSQ ( I ) ) 

DO  20  J=1.NS 
DMEVd,  J)  =0. 

20  EV(1.  J)=REAL(CEV(I,  J) ) 

DO  21  JJ=1.HS 
DO  22  1=1.  NS 
DO  S2  1  NS 

22  DMEVd,  JJ)  =DMEVd,  JJ)  ♦DDMd,  J)  »EV(J,  JJ) 
DO  24  1=1, HS 

24  DHORM(JJ)  =DHORM(JJ)-^EVd.  JJ)  MDMEVd,  JJ) 
21  CONTINUE 
DO  26  1=1, NS 

26  DNMd)  =DSQRT(1.  D  OO/DNORMd)) 

DO  28  1=1,  NS 
DO  28  J=1,NS 
28  Ed.  J)  =DNN(J)  «EVd,  J) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  BS,  B.  AND  BT  MATRICES  C 
C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  30  1  =  1,  NS 
DO  30  J=1,NS 

30  ET(J,  I)  =Ed,  J) 

CALL  GPROD(ET,  D.  BSS,  NS.NS,  K) 

DO  31  1=1, NS 
DO  31  J:1,K 

31  BSd,  J)  sBSSd,  J)/OMEGAd) 

DO  32  1=1,  N 

DO  32  J=1.E 
Bd,  J)  =BSd.  J) 

32  BT(J,  I)  =Bd,  J) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  C  INV.  AND  BCB  MATRICES  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CALL  AUGMNT (DK,  CCINV,  NS.  K,  IDGT,  LD) 

CALL  GPROD (CCINV,  BT,  CINVBT,  K,  K,  N) 

CALL  GPROD  (BS,  CINVBT.  BCB,  NS,  K,  N) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  IMSC  MODAL  GAIN  MATRIX  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  42  1=1,  N 
BETA (I )= BETA 1 /OMEGA (I ) 

42  ALPHA(I) =ALPHA1 
DO  43  1=1, N2 
DO  43  J=1,N2 

43  GM(I,  J)  =0. 
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DO  44  1=1,  H 
GH(2iil,  2iiI-l)=BETA(I) 

44  GH(2iiI,  2iil)  =ALPHA(I) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  TRANSITION  A  GAHHA  MATRICES  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  50  1=1, NS2 
DO  50  J=1,NS2 
GAHS(I,  J)  =0. 

50  PHIS(I,J)=0. 

DO  51  1=1,  N2 
DO  51  J=1,N2 
GAMH(I.  J)  =0. 

51  PHIH(I,J)=0. 

DO  52  1=1, NS 

PHIS(2iiI-l,  2«I-1)  :DCOS(OHEGA(I)  hDT) 

PHIS(2iiI-l,  2iil)  =DSIH(OHEGA(I)  kDT) 

PHIS(2iiI,  2x1-1)  : -PHIS  (2«I-1.  2«I) 

PHIS  (2x1,  2x1)  =PHIS  (2x1-1,  2x1-1) 

GAMS(2xI-l, 2xl-l)sDSIN(0MEGA(I) xDT)/OHEGA(I) 
GAMS(2xI-l,  2x1) : (1.  -DCOS (OHEOA( I ) xDT) ) /OMEGA ( I ) 
GAMS (2x1,  2x1-1) : -GAMS (2x1-1,  2x1) 

52  GAMS(2xI,  2x1) :GAHS(2xI-l,  2x1-1) 

DO  53  1=1, N2 

DO  53  J:1,N2 
PHIH(I,  J)  =PHIS(I,  J) 

53  GAMH(I,  J)  =GAMS(I,  J) 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  CLOSED  LOOP  EIGENVALUES  &  VECTORS  C 
C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

DO  60  1=1, N2 
DO  60  J=1,N2 
A(I,  J)  =0. 

ACLI(I,  J)=0. 

ACLN(I,  J)  =0. 

GHN(I,  J)  =.  0 
60  BCB2(I,  J)  =0. 

DO  62  1=1, N 
DO  62  J:1,N 

62  BCB2(2xI,  2kJ) =BCB(I,  J) 

CALL  GPR0D(BCB2,GH,  GHN,  H2,  N2,  N2) 

DO  64  1  =  1,  N 
A (2x1-1,  2x1) =OHEGA(I) 

64  A(2kI,  2x1-1)  =  -OHEGA(I) 

CALL  GHADD ( A,  GH,  ACL I ,  N2 ,  N2 ) 

CALL  GHADD (A,  GHN,  ACLN,  N2,  N2) 

CALL  EIGRF(ACLI,  N2,  LD,  IJOB,  EVALI,  EVECI,  LD,  WK,  lER) 
CALL  EIGRF(ACLN,  N2,  LD,  IJOB,  EVALN,  EVECN,  LD,  WK,  lER) 


•  fc 

»  ■  <  -  ^  ■  •  •  •  - 


o  o  o 


cccccccccccccccccccccccc 

c 

PRINT  OUTPUT  C 

C 

CCCCCCCCCCCCCCCCCCCCCCCC 
WRITE(6,  301)  H.  NS.  K 

301  FORMAT (//,  2X,  'N=*,  14,  4X,  'HS='.  14,  4X,  'E:',  14) 
WRITE(6,  302)  DT,  STIFF, DLEH,  RO, NUH 

302  FORMAT(//,  2X,  'DT=  " ,  E12.  4.  2X,  " STIFF:  ' ,  E12.  4,  2X, 
n'DLEN:  *,  E12.  4,  2X,  'R0:',E12.  4,  2X,  *NUH=',  13) 

WRITE (6.  303)  ALPHAl.BETAl 

303  FORMAT  (//,  2X,  'ALPHA:  ' ,  E12.  4,  4X,  '  BETA:  ',E12.4) 
WRITE(6,  304) 

304  FORMAT (//2X,  'D  MATRIX') 

DO  70  1:1,  MS 

70  WRITE(6,  100)  (D(I,  J) ,  J:  1,  K) 

WRITE (6,  305) 

305  FORMAT (//,  2X,  '  WSTINM  VECTOR  ') 

WRITE (6,  100)  (WSTINM(I),  1:1, N2) 

100  FORMAT(//,  2X,  7(2X,  D16.  7) ) 

WRITE(6,  306) 

306  FORMAT (//,  2X,  'STIFFNESS  MATRIX') 

DO  85  1:1, NS 

85  WRITE(6,  100)  (STIFM { I,  J) ,  J: 1 ,  NS) 

WRITE(6, 308) 

308  FORMAT (//,  2X,  'MASS  MATRIX') 

DO  86  1:1,  NS 

86  WRITE (6,  100)  (DMASS(I,  J) .  J: 1,  NS) 

WRITE(6,  311) 

311  FORMAT {//,  2X,  'C  INVERSE  MATRIX' ) 

DO  71  1:1,  K 

71  WRITE(6,  100)  (CCINVd,  J),  J:1,K) 

WRITE(6,  335) 

335  FORMAT (//,  2X,  'OPEN  LOOP  COMPLEX  EIGENVALUES') 

WRITE(6,  160)  (CW(I),  I:1,NS) 

WRITE(6.  336) 

336  FORMAT (//,  2X,  'OPEN  LOOP  COMPLEX  EIGENVECTORS') 
DO  72  1:1,  NS 

72  WRITE (6,  160)  (CEV { I,  J) ,  J: 1 ,  NS) 

WRITE(6, 337) 

337  FORMAT (//,  2X,  '  BS  MATRIX') 

DO  73  1:1,  NS 

73  WRITE(6,  100)  (BS (1,  J) , J= 1, K) 

WRITE(6, 313) 

313  FORMAT (//,  2X,  '  BCB  MATRIX') 

DO  74  1:1,  MS 

74  WRITE (6,  100)  (BCB ( I,  J) ,  J: 1 ,  N) 

WRITE(6,  314) 

314  FORMAT(//,  2X,  'E  MATRIX') 

DO  75  1=1, NS 

75  WRITE(6,  100)  (E( I,  J) ,  J= 1,  MS) 

WRITE(6,  315) 


V' V lu l IJ  H,  PJ IMUiJK.i'JiJ  JJ  J 
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315  FORMAT (//, 2X, ' IDEAL  OAIH  MATRIX') 

DO  76  1=1, H2 

76  WRITE (6,  100)  (QM( I,  J) , J= 1 ,  NE) 

WRITE(6, 317) 

317  FORMAT (//,  2X,  'HON- IDEAL  GAIN  MATRIX') 

DO  77  1  =  1,  H2 

77  WRITE(6,  100)  (GMH ( I,  J) ,  J= 1,  N2) 

WRITE(6, 316) 

316  FORMAT (//,  2X,  'COMPLEX  CLOSED  LOOP  EIGENVALUES') 
WRITE(6,  160)  (EVALN(I),  1=1,  N2) 

160  FORHAT(//,  2X,  3(2X.  2E16,  7)  ) 

WRITE(6, 320) 

320  FORMAT (//, 2X,  ' IDEAL  COMPLEX  CLOSED  LOOP  EIGENVALUES') 
WRITE(6,  160)  (EVALI (I),  1=1,  N2) 

WRITE (6.  324) 

324  FORMAT (//,  2X,  'COMPLEX  CLOSED  LOOP  MODAL  MATRIX') 

DO  78  1  =  1,  N2 

78  WRITE(6,  160)  (EVECN ( I ,  J) ,  J= 1 ,  H2) 

WRITE(6,  328) 

328  FORMAT (//,  2X,  'IDEAL  COMPLEX  CLOSED  LOOP  MODAL 
■MATRIX' ) 

DO  82  1=1, N2 

82  WRITE(6,  160)  (EVECI (I,  J) ,  J= 1,  N2) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C  C 
C  STRUCTURAL  RESPONSE  TO  MODAL  INPUTS  C 
C  (WSTS: STRUCTURE  MODAL  COORDINATES)  C 
C  C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

TIME=-DT 

KTIME=0 

147  KTIME=KTIME+1 
TIHE=TIME+DT 

IF  (KTIME.  EQ.  1)  GO  TO  510 
DO  200  1=1, N2 
WC(I)  =0. 

DO  200  J=1,N2 

200  WC ( I ) :WC ( I ) +GH ( I ,  J) nWSTS ( J) 

DO  202  1=1, NS2 
202  WCS(I)=0. 

DO  204  1  =  1,  NS 
WCS(2«I)  =0. 

DO  204  J=1,N 
CBW(J)  =0. 

204  WCS(2»I)  =WCS(2»I)  +BCB(I,  J)  «WC(2«iJ) 

DO  206  1  =  1,  NS2 
WSTS  ( I )  :  0. 

DO  206  J=1,NS2 

206  WSTS(I) =WSTS(I) ♦PHIS(I,  J) »WSTINS(J) 

DO  208  1=1, NS2 
DO  208  J=1,HS2 

208  WSTS ( I ) :WSTS ( I ) +GAHS ( I,  J) »WCS ( J) 


DO  210  1=1, NS2 
210  WSTINS(I) =WSTS(I) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c  c 

C  MODEL  RESPONSE  TO  MODAL  INPUTS  C 

C  (WSTM=MODEL  MODAL  COORDINATES)  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  250  1=1, H2 
WSTM(I)  =0. 

DO  250  J=1,H2 

250  WSTM(I)  :WSTM(I)  ♦PHIM(I,  J)  «WSTINM(J) 

DO  252  1=1, N2 
DO  252  J: 1,  N2 

252  WSTM(I)  :WSTM(I)  ♦GAMM(I,  J)  iiWCS(J) 

DO  254  1= 1, N2 
254  WSTINM(I) =WSTM(I) 

DO  256  1=1, r 
DO  256  J=1,H 

256  CBW(I)  =CBW(I)  +CINVBT(I,  J)  iiWC(2tiJ) 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  MODEL  &  STRUCTURE  C 

C  DISPLACEMENT  MATRIX  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

510  DO  508  J=1,NS 
DISP(KTIME,  J+1) =0. 

DISP(KTIME,  J+l+NS) =0. 

DO  507  1  =  1,  N 

507  DISP (KTIME, J+ 1 ) =  DISP (KTIME,  J+ 1 ) +E ( J,  I ) «WSTM ( 2 « I - 1 ) 
DO  509  1=1, NS 

509  DISP (KTIME,  J+l+NS) =DISP (KTIME,  J+l+NS) +E(J,  I) « 
«iWSTS(2«I-l) 

DISP (KTIME, J+l+NS) =40.  «DISP (KTIME,  J+NS+1 ) 

DISP  (KTIME,  J+ 1 )  =  40.  »DISP  (KTIME,  J+ 1 ) 

508  CONTINUE 

DO  511  J=1,K 

511  DISP (KTIME,  J+l+NS+NS) =40»CBW(J) 

DISP (KTIME,  1) :TIME 

IF(TIME.  LT.  10.  OD  00)  GO  TO  147 

CALL  MPLOTdO.  D  00,  2.  OD  00,  KTIME,  NS,  NUM) 

STOP 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  MATRIX  AUGMENTATION  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  AUGMNT(A,  B,  N,  M,  IDGT,  LD) 

IMPLICIT  REALi*8(A-H,  O-Z) 

DIMENSION  A(40, 40) , B(40, 40) , A1 1 (40, 40) , A12 (40, 40) , 


98 


kA2 1 (40.  40) .  A22 (40,  40) .  AZ (40.  40) ,  WEAREA ( 500) . 
IIAI22  (40.  40) .  AIA(40.  40) 

NB=N-H 

IF  (HB.  EQ.  0)  THEN 
DO  5  1:1, H 
DO  5  J:1,H 
5  B(I,  J)  :A(I,  J) 

EHDIF 

IF  (HB.  EQ.  0)  GO  TO  10 
DO  1  1:1, H 
DO  1  J:1,H 

1  All  (I,  J)  :A(I.  J) 

DO  2  1:1,  H 

DO  2  J:1,HB 

2  A12(I,  J):A(I,H+J) 

DO  3  I:1.HB 

DO  3  J:1,H 

3  A21 (I,  J) :-A(H*I,  J) 

DO  4  I:1.NB 

DO  4  J:1.NB 

4  A22(I,  J)  :A(H+I,  H+J) 

CALL  LIHV2F(A22,  HB.  LD,  AI22,  IDGT,  WZAREA,  lER) 
CALL  GPROD(AI22.  A21.AIA.  HB,  HB.  H) 

CALL  GPR0D(A12,  AIA.AK.M,  NB.  H) 

CALL  GMADD(A11,  AK.  B.H,  H) 

10  CONTINUE 
RETURN 
END 

CCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  MATRIX  PRODUCT  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  GPROD (A,  B,  C,  M,  N,  K) 

IMPLICIT  REAL«i8(A-H,0-Z) 

DIMENSION  A (40, 40) .  B(40.  40) ,  C(40.  40) 

DO  1  I:1,M 
DO  1  J:1,K 
C(I,  J)=0. 

DO  1  L:1,N 

1  C(I,  J):C(I,  J)+A(I,L)»B(L.  J) 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  MATRIX  ADDITION  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  GHADD(A,  B,  C,  N.  M) 

IMPLICIT  REAL«8(A-H,0-Z) 

DIMENSION  A (40,  40) ,  B(40,  40) ,  C(40,  40) 

DO  1  1:1, N 
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I  C(I,  J)=A(I,  J)+B(I,  J) 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  PLOTTING  ROUTINE  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  HPLOT(XHAa, XT,  H.  EH,  EPT) 
IMPLICIT  REALii8(A-H,0-Z) 
COHHON/COHCIZ/A ( 5002,  46 ) ,  YH ( 1 0 ) ,  YT ( 1 0 ) 
XSCALE=10.  /lO. 


XORG=i. 


Y0RG=5.  5 

CALL  PLOTS  (0,0,0) 

DO  5  1  =  1,3 

CALL  PLOT(l.  0,  0.  ,  -3) 

CALL  PLOT(0.  ,  11.  ,  2) 

CALL  AX ( XHAG,  XT.  YH ( I ) ,  YT ( I ) ,  XORG,  YORG ) 
YSCALE=3.  5/YH(I) 

K:3 


IFd.EQ.  1)  NrEPT+l 
IFd.EQ.  2)  NsEH+KPT+1 
IFCI.EQ.  3)  NsEH+KH+KPT+l 
DO  1  J=1,M 

4  X=A(J.  1) kXSCALE+XORG 
Y=A( J,  N) «YSCALE+YORG 

3  CALL  PLOT(X,  Y,K) 

IF(K.  EQ.  3)  Kr2 
1  CONTINUE 

6  CALL  PLOTdO.  5,  0.  , -3) 

5  CONTINUE 

CALL  PLOT(0.  ,  O.  ,  999) 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C  C 

C  PLOTTING  SUBROUTINE  C 

C  C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  AX (XMAG,  XT, YHAG,  YT.  XORG,  YORG) 
IMPLICIT  REALi«8(A-H,0-Z) 

XH=10. 


YH=  3.  5 

CALL  PLOT (XORG,  YORG- YH-.  5,  3) 
CALL  PLOT (XORG,  YORG+YM+.  5,  2) 
Y=YH«YT/YHAG 
YY=YM 


YI=YH+. 1 

1  CALL  PLOT(XORG,  YORG+YY, 3) 
CALL  PLOT (XORG-.  1,  YORG+YY,  2) 


IF(YY.  GE.  -YI)  GO  TO  1 

X:XHiiXT/XHAG 

XX=XH 

CALL  PLOT(XORG,  YORG,  3) 

CALL  PLOT(XORG+XH+.  5, YORG,  2) 
2  CALL  PLOT(XORG+XX.  YORG.  3) 
CALL  PLOT(XORG+XX,  YORG-.  1, 2) 

xx=xx-x 

IF (XX.  GT.  .  1)  GO  TO  2 

RETURN 

END 

/» 

//GO.  DATA  DD  « 

(  Place  Data  Here  ) 

/« 

//GO.  PLOTPARH  DD  « 

&PLOT  XHAX=120  &EHD 
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